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Abstract In this present work we consider a falling body through the air 
under the influence of gravity. In particular, we consider the experimental 
data based upon the free fall of six men in the atmosphere of the earth. 
In order to describe this process we employ a nonlocal dissipative force. We 
show that our description, by using an exponential memory kernel, can flt 
the experimental data as well as that obtained by a local dissipative force. 
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Recently, the idea of employing a nonlocal dissipative force has captured so 
much interest in the scientific community. For instance, in diffusion processes, 
generalized Langevin equations have been employed to describe anomalous 
diffusion regimes, including both subdiffusion and superdiffusion. According 
to Kubo [1] the frictional memory kernel is related to the form of the corre- 
lation function for random force when the system is in equilibrium state. As 
can be noted, anomalous diffusion transport is present in diverse physical sys- 
tems, e.g., it can be found in porous media [2l|3j, polymers [11[5], amorphous 
semiconductors [6l[7] and composite heterogeneous films [8]. 

Approaches based on the fractional derivatives have also been employed 
to describe, for instance, anomalous diffusion transport [^llOj and viscoelastic 
systems [TT | fT2 l [T3| . The fractional derivatives are also nonlocal and they can 
be related to memory effects. Recent works on the theme, have revealed that 
a fractional derivative can behave as a damping term [9lll0j. In this respect, a 
simple fractional oscillator has also been employed to investigate the physical 
behavior of a fractional derivative, and it supports this point of view [14pi5j. 

In this work we consider a simple process involving a falling body through 
the air. We also consider an observed process based upon the free fall of six 
men in the atmosphere of the earth. As it is well-known the usual description 
employs the Newton second law with a dissipative term in function of the 
velocity of the falling body. For this observed process the usual model, with 
the dissipative force proportional to w^, is employed. This phenomenological 
model is nonlinear. However, other descriptions can also be used to describe 
the process and they can fit the observed data equally well. In our description, 
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we work on the assumption that the dissipative force is nonfocal which is 
related to the memory effect. We analyze our model by employing short and 
long memory kernels. Our results are compared with the result of the usual 
approach and experimental data. 



2 Local and nonlocal models 

When an object falls from rest through the air surrounding the earth, it 
experiences a resisting force that opposes to the relative motion in which 
the object moves relatively to the air. Experimentally, this resisting force is 
related to the relative speed v. For slow speeds the resisting force is in magni- 
tude proportional to the speed. This last assumption is usually employed to 
describe diffusion processes such as in the Langevin equation [16j . However, 
in other cases it may be proportional to the square (or some other power) of 
the speed. Applying the Newton second law pTifTS] we haverwe have: 

— mg — mF{v) , (1) 

where g is the free-fall acceleration and F{v) is a function of the velocity. 
We note that for F{v) ^ v'^, Eq. (1) is nonlinear and it corresponds to the 
well-known Riccati equation [19j . In particular, for F{v) — Kiv/m, we have 
the following solutions for z and v: 



v{t) 



mg 



ci exp 



(2) 



and 
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z{t)=ZQ- — + —{^gt + ciexp^-—t^j, (3) 

where ci = mg/Ki—vo, zq is the initial position and t;o is the initial velocity. 

For F(v) = Kiv^ jm, we obtain 



1 - C2 exp [— yicj 

and 



^, /(I -C2CXP -pi ) \ 

p V(l-C2)2exp[-pt] J 
where b"^ = mg/K2, p = 2bK2/m and C2={vo — b)/(vo + b) . 

For these cases, Eq. (1) shows that the resisting force of an object increases 
with the increase of the relative speed v. Consequently, when the body falls 
enough, it reaches a constant terminal velocity. For solutions (2) and (4) the 
terminal velocities are obtained by taking t — > oo, and the results are: 



"-S? (e) 

and 



Now one considers the above system by assuming a nonlocal dissipative 
force given by 

F{v)= I '){t-T)v{T)&r , (8) 
Jo 
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where j{t) is the memory kernel. We note that Eq. (1) is Hnear with F{v) 
given by (8), and its solutions can be obtained by using the Laplace trans- 
form. For 7(f) ~ S{t) one has a local dissipative force with F{v) ~ v. For a 
nonlocal memory kernel we choose the exponential and power-law functions. 
We see that these functions have short and long tails, respectively. Our idea, 
with these functions, is to analyze the influence of memory function on the 
system described by Eq. (1) for short and long memory effects. 

For exponential memory kernel 7(f) = 706"*^* we obtain the following 
solution: 



qA 

z = zo+A + —t■ 
lo 



A 



70 



Vo 



70 2 V " 70 2 J 



g\ A 
70 



e ^% 70 = 



9A 

z = zo + A+ —t+ 
70 



[Vo- — - -A) — ^—^ - A cos(Bit) 
V 70 2 / Bi 



e ^\ 70 > 



4 ' 



(10) 



(11) 



(t) = — + (vo - —) cos(Bit)t 
70 V 70^ 

A / gX X 



ABi (vo-- 

2Si V ° 70 2 



A 



sm{Bit)t 



qX 

z = zo + A+ —t+ 

70 



/ 9X X \ sinh(i?2i) , , /^^ \ 

[vo- — - -A) A L- A cosh 

V 70 2 / B2 



7o>y , 



-At A2 

e ^% 70 < 



(12) 



(13) 



and 
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AB2 + —\Vo 

2B2 V 70 



v(t) = ^ + (vo~ —) cosh(B2t)e~^ 



gX 
70 

gX X 



smh{B2t)e 2*, 70 < 



4 ' 



(14) 



where A= {g + Xvo- 5^^/70) /70, = Vlo - A2/4 and B2 = y^A2/4-7o. 

For the three cases we can verify that v{t = 0) = vq and they give the 
same expression for the terminal velocity which is given by 



70 



(15) 



Now, it is more convenient to write Eqs. (9), (11) and (13) in terms of the 
terminal velocity 



voV 31/2 

z = zo + — — + Vt 

9 4 51 



4 g 



z = zo+A+Vt- 
and 



+ 2 ) 

(.o-^-^a)!^-Acos(B.O 



e V ' 



70 



A2 



, (16) 



-4* 



e"^% 70 > (17) 



z = zo + A + Vt+ 



X \ smh(B2t) 
vo-V- -A) - ^cosh 



e 2% 70 < 



A^ 



(18) 



where A = V/X+voV/g-V^ /g, Bi = ./gX/V-X^/A and B2 = v'AV4 - 5A/V. 
We note that 70 does not appear explicitly in these last expressions. In par- 
ticular, the solution for 70 = ^ only contains one unknown parameter such 
as in the usual approach. For two other cases the free fall can be described 
by one unknown parameter A which can be determined by the experimental 
data. 

In order to obtain some insights on the solutions of the nonlocal approach 
(9)-(14) and the solution of the usual approach given by Eq. (5) let us fit the 
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experimental data from the free fall of six men in the atmosphere of the earth, 
from a maximum altitude of 31,400 to 2,100 feet which was accomplished in 
116 seconds [19) . The average weight of the men and their equipment was 
261.2 pounds and the average velocity of the lowest point corresponding 
to a drop of 29,300 feet in 116 seconds is V ~ 251.4 feet per second. The 
experimental data are shown in Table 1 with the values of H — z, H = 31, 400. 

To find the velocity of the lowest point V we set vq = and zq — 0. We 
obtain from Eq. (5) the following approximate expression 



Substituting z ^ 29,300, t = 116 and g ^ 32.2 wc find V2 = 265.7 ft/s. 
This value is somewhat larger than the experimental result. In the case of 
F{v) V, the terminal velocity is larger than V2, thus this model should be 
put aside. For the solution (16) we find the following approximate expression 

3y2 

z = vt--— . (20) 
4 <? 

Note that the expressions (19) and (20) are similar and the values In 2 and 
3/4 are close. From Eq. (20) we find V = 266.89 ft/s. The parameter A can 
be estimated from (15) and we find A = 0.483 s^^. 

In Fig. 1 we show the comparison of the observed values with those ob- 
tained by the solutions (5) and (16). It is clear that the solution (16) gives 
the adjustment as well as that of the usual approach. In Fig. 2 we show the 
behavior of the solutions (5), (16), (17) and (18). For the solution (5) the 
terminal velocity is V2 — 265.7 ft/s, while for other solutions the terminal 
velocity is V = 266.89 ft/s. The plots show that the curves are close to each 




(19) 
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other. Thus, aU of them can be used to describe the free fall system equally 
well. 

Next, we consider a long tail memory kernel given by j(t) = 7it~" for 
< a < 1. The solutions for z{t) and v{t) are given by 



z{t) = zo + gt^E2-c.,3{-l*t^'") + votE2-aA-l*t^~'') (21) 

and 



v{t) - gtE2-cA-l*t ) + voE2-c.A~l*t ) > (22) 

where 7* — 71/^(1 — a) , r{x) is the gamma function and £"^,1/ is the general- 
ized Mittag-LefHer function defined by Ef,,^(z) = Y.'^^oi^T / r{nn + [20l 
[21] . For i = we verify that v{t = 0) = vq. In order to verify the terminal 
velocity we should analyze the behavior of Eq.(22) for large time. To do so, 
we use the following asymptotic approximation: 



E,A^) jTT^ r . (23) 



Then, we obtain from (22), for vq = 0, the following approximation: 

Ht) - 4-^ . (24) 

7*i [aj 

This last result shows that Eq. (22) can not give a terminal velocity. Then, a 
long time memory such as given by the power-law memory kernel, based on 
Eqs. (1) and (8), is not appropriate to describe the free fall system through 
the air. 
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3 Conclusion 



In summary, we have investigated the faUing body problem by using the lo- 
cal and nonlocal dissipative forces. We note that nonlocal dissipative forces 
are frequently used to describe diffusion processes by a generalized Langevin 
equation. Then, this kind of approach represents a well-known generalization 
of the usual approach. We have shown that the results of a nonlocal dissipa- 
tive force given by the exponential memory kernel can fit the experimental 
data as well as the result obtained by the local dissipative force F{v) . 
In contrast with the nonlinear force F(v) ^ w^, the nonlocal dissipative force 
(8) gives a linear differential equation. On the other hand, the power-law 
memory kernel does not give a terminal velocity and consequently it alone 
can not be used to describe the free fall system. These results suggest that 
a short memory effect in the nonlocal approach may be appropriate to de- 
scribe the free fall system. The presence of memory effect in the system may 
be related to the fact that when the body moves through the air, it displaces 
the particles in its immediate vicinity. Then, the surrounding flow field is 
altered and acts back on the body, resulting in an action which depends on 
the past motion of the body. This feedback process which leads to a memory 
effect can be reflected to the motion of a particle in a fluid [52]. Further, we 
have also verified Eq. (1) with an additional term for the dissipative force 
given by F{v) — Kiv/m + 70 exp(— A(t — T))t;(r)dr. For this case we can 
improve the adjustment, i.e., we can fit the experimental data as well as the 
result of the usual approach with a lower terminal velocity {V — 265 ft/s). 
In Fig. 3 we show the comparison of the result with that described by the 
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usual approach. In the case of a fractional derivative, it can also be employed 
to replace the classical Newtonian force, but it does not lead to improve the 
result obtained by the usual approach significantly [23j . 
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Figure Captions 



Fig.l - Adjustment of the experimental data with the expressions (5) 
(dotted line) and (16) (solid line). The terminal velocity for the expression (5) 
is given by V2 = 265.7 ft/s, whereas the terminal velocity for the expression 
(16) is given hy V = 266.89 ft/s. 

Fig.2 - Plots of the expressions (5), (16), (17) and (18). The terminal ve- 
locity for the expression (5) is given by V2 = 265.7 ft/s, whereas the terminal 
velocity for the three other expressions is given by ^ = 266.89 ft/s. 

Fig. 3 - Adjustment of the experimental data with the expression (5) 
(dotted line) and the result obtained by Eq. (1) with the dissipative force 
given by F{v) = K\v/m + 70 /g exp(— A(t — T))v{T)dT (solid line). The ter- 
minal velocity for the expression (5) is given by V2 = 265.7 ft/s, whereas 
the terminal velocity for the second case is given hy V = 265 ft/s, where 
Ai = A -I- Ki/m and 70 + XKi/m = Af /4. 
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31,400 


46.7 


20,550 


88.1 


9,400 


9.9 


30,780 


51.3 


19,400 


92.7 


8,140 


14.5 


30,200 


55.9 


18,100 


97.3 


7,080 


19.1 


28,850 


60.5 


16,600 


101.9 


5,700 


23.7 


27,700 


65.1 


15,150 


106.5 


4,450 


28.3 


26,150 


69.7 


14,070 


111.1 


3,170 


32.9 


24,600 


74.3 


12,800 


116.0 


2,100 


37.5 


23,200 


78.9 


11,650 






42.1 


21,750 


83.5 


10,400 







Table 1 Estimated values of H — z for the fall from an altitude of 31,400 to 
2,100 feet, with z = Q when t = Q. 



